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Abstract 

Motivated by the recent experimental data about candidates for glueball from different processes, 
we discuss in this paper the production of ++ glueball from double diffractive scattering at momentum 
transfer |f|<^lGeV 2 in high energy p + p(p) collision. We apply the phenomenology of Pomeron (iP) of 
Donnachie-Landshoff, the field theory model of IP of Landshoff-Nachtmann and the relevant calculating 
approachs. We assume while IP coupling with glueball, the ++ glueball can be considered as a bound 
state of two non-perturbative massive gluons. We evaluate the dependence of cross section for ++ 
glueball production on system energy y/s and show that it could be tested experimentally. 



PACS number(s): 12.40Nn, 12.39.Mk, 13.85.Lg 



1 Introduction 



Since the color SU (3) gauge theory possess of self-coupling characteristic between gluon fields, it had pre- 
dicted long before the quantum chromodynamics (QCD) became the basic dynamics of strong interaction, 
that there would exist states in hadronic spectroscope which are formed from purely gluon field — glueball |l|]. 
Though now there are still lack of strict theoretical proving and conclusive experimental evidence to confirm 
its existence, in recent couple of years there are a lot of reports from relevant data analysis that several 
candidates of glueball and/or states mixing with quark pair have been observed experimentally It is well 
believed that glueballs should easily be produced in processes with rich source of gluon constituent, thus 
experimentalists concentrate their attention to seek glueballs on pp annihilation ||, J/ip radiative decay Q], 
and high energy central h—h collisions Glueball production from J /if) radiative decay has been analyzed 
in a model independent manner by Close et. al. || , but yet there has not been any analysis in central h—h 
collision and in pp annihilation. 

At high energy central h — h collisions the increase in cross section with increased center of mass energy 
JtJ is consistent with the double pomeron exchange (DPE) machanism , which was predicted to be a source 
of gluonic states A large glueball production cross section in the central region is predicted by Gershtein 
and Logunov, they related the rise of total cross section with increasing energy to the exchange of glueballs 
in the t— channel or to the collising of the soft gluon seas of the interacting particles [||. These observations 
tell us that glueball production in high energy h - h central collisions would be intimately connected with 
the non-perturbative mechanism . 

In this paper, we discuss the glueball production from double diffractive scattering in high energy p-\-p(p) 
collisions, 

p + p(p)^P+p(p) + G, (1) 

as a first attempt, we consider J PC = ++ glueball only. 

Our approach is on the basis of the phenomenology of IP from Donnachie-Landshoff ]l(J , the field theory 
model of IP from Landshoff-Nachtmann|ll|. Based on the above approach, a more refined and sophisticated 
model of diffractive scattering had proposed by Cudell and Hernandez [ ]l2| , concentrating on Higgs diffractive 
production in pp collision p + p — > p + p + H, they obtained a result close to that of Bialas and Landshoff|l3j 
using the L-N model of DPE. 

Since the glueballs are produced through DPE, it could believed that they are laid predominately in the 
central region of final state rapidity distribution and with large symmerty rapidity gaps from each of the 
final diffractive protons ( or antiprotons) which are the characteristic properties of DPE processed as Eq.(l) 
. In principle, it should be easily detected and according to it we could distinguish it from all other strong 
backgrounds. We will come back to this point in last section. 

Of how to link L-N field theory of IP with the phenomenology of Regge behaviour, we follow the approach 
of Bialas and Landshoff jl3| . Since the information of coulping of ++ glueball with nonperturbative gluon 
coupling is very scarely , we adopt the color singlet approximation fll4|] , which always used in heavy qurak 
physics. Under this approximation the coupling structure is Bg a ^S e h , where a e etc. are the Lorentz and 
color indexes of the two non-perturbative gluons respectively , the parameter B can be estimated from the 
branching ratio of J/ip radiative decay to ++ glueball. 

In the following we describe in detail our calculating scheme in Sec. II, we give the formalism and derive 
the coupling constant B of the two nonperturbative gluons with ++ glueball in Sec. Ill, our results and some 
related discussions are given in last Sec. IV. 



2 Calculating Scheme 

In phenomenology of IP given by Donnachie-Landshoff |l0| , the nucleon is to be treated as a bound state 
composed of three clusters, they are formed from non-perturbative QCD effects and their size are much 
smaller than nucleon radius. The core of each cluster is a valence quark, with sea quarks and gluons 
surrounded. Furthermore D-L have demonstrated that in coupling with nucleon, to a good approximation, the 
IP behaves like a isoscalar ( C= + I) photon. So it means that the coupling of IP with a nucleon is actually 
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that of IP with the three clusters (constitude quark ) in an incoherent manner. In this picture IP model one 
does not need to consider the interaction of IP with any other parton (sea quark or gluons) in the nucleon. 
With this model Donnachie and Landshoff have successfully explained a lot of experimental effects in high 
energy soft strong interaction processes pl| . 

In discussing glueball production from DPE , we adopt the same point of view: we should not start such 
a problem from parton model and perturbative QCD, but instead from non-perturbative QCD at all and 
assume the glueball to be a bound state of two small flavorless clusters (or "constituent" gluons) which are 
formed from gluons by non-perturbative QCD effects, and are massive. Like the nucleon case, we assume 
the coupling of IP with glueball is actually that of IP with these "constituent" gluons. Also, we do not need 
to consider the interaction of IP with any other massless gluons in the glueball. 

From the status both of the available data for elastic diffractive scattering and the effectuality of D-L IP 
model in high energy h-h collisions, we restricted in this paper the momantum transfer of scattered proton 
(antiproton) to to be small ,e.g \ti\ <l.0GeV 2 . The whole picture is sketched in Fig.l. 

From Regge pole theory |l6[] and IP model of D-L, when s < si, S2 < m^r, and \ti\<lGeV 2 , the 
asymptotic form of the amplitude for Fig.l is 

(±)«-(t I )-l(±)^(t2)-l F (G) (aip(tl))aip(t2) ^ )i?l(tl)i , l(t2)7A87 A (2) 
S2 Si 

where tj = (pi — p^) 2 , z = f^r, and m is the ++ glueball mass, aip(t) is Regge pole trajectory of IP. From 
data analysis when \t\ is small, opp(i) w a/p(0) + a' lp t = 1.086 + 0.25GeV~ 2 t Pi(t) is the Dirac form 
factor of proton. The direct product of gamma matrix shows that the external proton lines will be traced 
when calculate the cross section. 

Since the diffractive scattering condition requires 

(3) 

s s 

and usually assume S — 0.1, hence in the asymptotic form Eq.(2) one could neglect lower powers of (^-) and 
(-^)- F' G ) (o!/p(ti), aip{t 2 ), z) is the IP — IP — G vertex function, it general structure is known ]l7t, but its 
concrete form must be further considered. 

Since IP-nucleon coupling should be considered as pomeron incoherently coupled with three constituent 
quarks, and according to field theory model, the IP exchange is just the exchange of two non-perturbative 
gluons [p"l| Jill , thus Fig.l should be concreted as Fig. 2. 

In the field theory model of pomeron, the IP — IP — G vertex function F^ G \aip{ti), otipfa), z) for 
\ti\ <1.0GeV 2 can be calculated at by the sum of Feynman amplitudes of Fig. 2 which as shown in Jl3j is 
equal to its s— channel discontinuity of the first diagram, thus as Fig. 3. This simplify our calculation greatly. 
The black blob in Fig. 3 is the vertex function of ++ glueball with gluons of IP which will be discussed in 
next section. 

In the following we first get the forward diffractive scattering amplitude for Eq.(l) in L-N field theory 
model of IP, then identify it with the corresponding Regge behaviour formula Eq.(2) at ti = to obtain the 
normalized value of IP — IP— glueball vertex function f ( G ) (app(ii), ctipitz), z)\ tl=t2=0 . 1 After this has 
done we continue Eq.(2) to \ti\ ^ and discuss DPE production of glueball. 



3 Formalism and Coupling Constant B 

Using the Sudakov variables to calculate the loop integral of Fig. 3 is very convenient, here 



xpi | VP2 
s s 



1 Since in D-L model, ctjp(t)\t=o = 1 + £ ~ 1.086, so it seems having some inconsistent with this identification, which has 
explained in jlij. Furthermore, this problem has studied by Ross [|3 who using a hybrid model in which one uses PQCD to 
treat interactions between gluons, but takes a D-L type non-perturbative propagator for the gluon. In the leading logarithm 
approximation, by summing a dominant subset of diagrams, the phenomenolog required Regge behaviour is obtained. 
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p 1 = xipi h h vi 

s 

P2 = — + V2P2 + V 2 (4) 

We have put the light quark mass equal zero, the v, v\, and vi are transverse to p\ and p 2 and so effectually 
are two dimensional, v 2 — — v 2 ,t>| = — v 2 . Then 



h 


= (pi 


-ft) 2 = 




t 2 


= (P2 


"^) 2 = 


— vi/ifa 


Sl 


~ 8(1 


- Jfc) 




S2 


~ *(1 


- Xl) 





(5) 

As mentioned above, we first need the s— channel discontinuity of Fig. 3 at t\ ~ t% ~ 0, so we set 
v\ f=a V2 ~ 0, then 

d 4 M( gi )<%i) ~ ±-Jdxdyd 2 v6(y-v 2 )5(-x-v 2 ), 
J d*p>M 2 ) ~ i J dx 1 dy 1 S(x 1 yi)d 2 v 1 , 

J d^p^Sip^ 2 ) ~ \ J dx 2 dy 2 8(x2y2)d 2 v 2 . (6) 

Here we explain how we can get from Fig. 3 the direct product from 7a ® 7 A as in Eq.(2). The upper 
quark line in Fig. 3 has gamma matrices 7^7 • 91 7a, for large s its asymptotic form equivalent to 

2<7\ (7) 

because in calculating the differential cross section, we also need to multiply 7 • p\ on the left and 7 • p[ on 
the right of this expression. When these are included, the difference between (7) and the original expressions 
gives contribution to cross section are of order 6. Similarly, from the lower quark line we obtain 7a- 
Now from Fig. 3, the amplitude of q — q(q) diffractive scattering through DPE at tx = t% = is 

M (G) 7 A ®7 A , 

M<°> =M- [ d 2 vW^D(-v 2 )D{-x 1 x 2 )D{-y2V 2 ), (8) 
2-k z s J 

where D(q 2 ) is non-perturbative gluon propagator of IP, g is the coupling constant between non- 
perturbative gluon with quark, 

g 2 D(q 2 ) = -A exp(^), (9) 
A* 

where /i = l.lGeV, A 2 p 2 = 72-K0 2 ,(3 2 = 3.93GeV~ 2 is IP— nucleon coupling constant in the pomeron D-L 
model @. 

Using the approximation matrix form Eq.(7), we found from Fig.3, the W^ G) in Eq.(8) , 

= q£A$(g, (10) 

where vertex function A^J connecting the gluons of IP with ++ glueball. The simplest diagrams for this 
vertex function shown in Fig. 4 and Fig. 5. 
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In deriving these vertex functions, for calculating simplicity but not necessary we have used color singlet 
model approximation (l4[ for the vertex, which require gluon lines cutting by vertical dashed line in these 
figures are limited on mass shell. Thus from Fig. 4 we get 



A %c = 2BF^ aX (k 1 , q ,~^P)D(q 2 )g^F^ p (-^k 2 ,~q)gP x -f be J hcg -6 eh , (11) 

where 

F„ l ^{ri,r 2 ,r 3 ) = {r l - r 2 ) fi3 g ttltt2 + (r 2 - r 3 ) M1 e^ 2AI3 + (r 3 - r x )^g^ 

ri+r 2 +r 3 = 0. (12) 



From Fig. 5 we get 



From Eqs.(10) to (13) we get 



\v,bc - -^Bg^fbegfhcgSeh- (13) 



W G 3 = yq ' . (5 TO V + IOtoVv 2 + 16m 2 sv 4 + 15s 3 v 2 a;i + 15s 3 v 2 y 2 - 10s 3 v 2 

+6s 2 v 4 xi + 6s 2 v 4 y 2 + 4s 2 v 4 + 12sv 6 xx + 12sv 6 y 2 + 24v 8 ) 



q 2 = — (m 2 s + 2sv 2 x 1 + 2sv 2 y 2 + 4v 4 ) 
4s 

1 v 4 
W Gi = -12£(~s + — ) 

2 s 

W (G) = W G 3 + W G± ( U j 

Form Eqs.(8)(9)(14) , we get M^ G \ for Fig.3, which is F^ G \aip{t-s), ajp(t 2 ), z)\t 1= t 2 =o as we explained 
in section II . 

For ti, t 2 0, the general structure of F(aip(ti), aip(t 2 ), z) can be analyzed from Regge theory and this 
has been done by Drummon et. a/. |l7|] . They have shown that when \t 2 \ are small, F(aip(ti), aip(t 2 ), z) 
are finite for any z. This point is very important for us since the D-L IP model, in strictly speaking, is valid 
only for \t\\, \t 2 \<lGeV 2 . But we can also see from Eqs.(2) and (3), all contributions from large |ti|, \t 2 \ 
are higher powers of 8, and they can be neglected safely. So, for ti, t 2 ^ 0,we have 

-F(a/p(ii),a/p(i2),z)|^|,|«2|<iG e v 2 ~ F(ajp(ti), aip(t 2 ), z)\ tl =t 2 =o = M {G) (15) 

Thus the amplitude T' G ) for glueball production in p — p(p) high energy double diffractive scattering from 
D-L and L-N IP model and formalism, in a good approximation, can be got from Eqs.(2) and (15) and is 

T iG ) = 9Af( G )(-) t " ff, ( fl )- 1 (-) aff> ( i2 )- 1 F 1 (t 1 )F 1 (i 2 )7 A <g> 7 A . (16) 

where factor 9 comes from 3 quarks in nucleon, F\(t) is the Dirac form factor of proton, using the exponential 
approximation for at small \t\ : F\(t) ~ e xt , A ~ 2GeV 2 , the cross section formula is 



. J^f / -^1 |^((l - - y 2 ) - =L) 

2(27r)° m z J x\ y 2 s 

• /d 2 Wl d 2 W 2 (l-Xi) 2Q ' V i/^(l-y2) 2Q ' vi/a2 . e ^Kvl/xi+vl/y 2 ) (17) 



where F G = | is the color factor. 
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Let us now consider the coupling constant B of constituent gluons with ++ glueball, which could be 
fixed by radiative decay mode J/ip — ► 7 + /o(1500) as shown in Fig. 6, ^ where 

Pi = t;{pj-Pg) 

Pj = k- ^pj 

Pk = ^(PG-Pj) 

.91 = 92 = ^Pg (18) 

It is well know that the J/ip is dominatly an S- wave state and the non - relativistic quark potential 
mode is very successful in describing the static properties of the J ftp. Inclusive quarkonium decays into light 
hadrons are accessible to PQCD[|22| . For heavy- quarkonium, the annihilation time ~ is much smaller 
than the time scales relevent to QQ binding. Hence the short - distance annihilation of the QQ pair can be 
separated from the long- distance effect of QQ binding. In our J /if) — > 7 + G case, the long- distance effects 
of QQ and gg binding can be separated, the gg binding effect is included in the non-perturbative coupling 
constant B of constituent gluons with the ++ glueball . 

Using the color singlet model approximation (l4| in J/ip vertex as usual, we get the amplitude of J/ip — > 

7 + /0 M J = 2A e q g%B(M$ p + M<& + M^g^e^ (7) (19) 

where e M (fc) is the polarization vector of photon, e q = |e is the charge of charm quark, coefficient Ao fixed 
from Tj^ e+e -: 

Ale q = r ^-^ ( 20 ) 

and 

1 7 • pj + m c 7 • + m c 
o7 ' PJ + m c 7 M — 2 — > — 2 2~ Tp 

2 Pj ~ m c Pi — m i 

A \ l-Pk+m c -7 -pj +m c 
b7 • Pj + m c )^ v — 2 — 7 P — 2 — 7/* 

2 Pfc — "^c — m c 

/I . x J-Pk+m c j- Pl + m c 

2 P k - m i Pi~ m c 

Put Eqs.(18) and (19) into (17), after sum and average with initial states and final states respectively, 
we get 

^ 7 n 2 „ n o « 1024(m 6 - to 2 to 4 + \\m 2 mS + m 6 ,) 

53 |M''| 2 = -A^B e gQ ^ ^ ^ J -±, (22) 

3 * rrij{rnj — m z y 

then decay width of J/?/' — > 7 + /o is 



M (a) = 


—=Tr[j 
V2 


M< 6 > = 

flVp 


— Tr[7 
V2 


M (c) = 


— =Tr[7 
V2 



'.^ + /o = ^^^El A/J | 2 ' ( 23 ) 
lbnrrij z — ' 



where color factor F/ = |. 



Put Eqs.(18), (20), and (21) together, we get 

B 2 = r J^7+/o 9nam e j(m 2 } - to 2 ) 

rj^ e+e - 128^(to^ + 11to 4 7 to 2 -to 2 7 to 4 + to 6 )' 1 j 

Decay width Fj^ e + e - and r,/^ 7+/o can be found in p3]|, ^ J ^ T + /o = g^^.j so as go and m are fixed, one 

gets coupling constant B. For f| = 50 ; n [mj/A] ' A = 0-20GeU,TO = 1.5GeV , we getB 2 = 3.8 x lO^GeU 2 . 

2 Close argue in |2l|] that the scalar /p(1500) may be a glueball— qq mixture, here for simplicity we assume it is a pure 
glueball. 
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4 Results and Discussions 



We first give a brief comments on the parameters used in our paper. 

1. ++ glueball mass m . 

From several models {e.g. lattice QCD, bag model, potential model, and QCD sum rule) analysis, 
the lowest state mass of ++ all are fixed at 1.5 ~ 1.7Gey[|TJ. Experimentally, relevant data also 
manifest a clear signal of ++ resonance state at about l.hGeV in the central region of high energy 
p — p collision, thus we let m = l.5GeV. 

2. Nonperturbative coupling constant g. 

In the Abelian gluon theory of Landshoff and Nachtmann, the nonperturbative gluon only couples to 
quark, thus g 2 D(q 2 ) always appear together and can be normalised by the constant A in Eq.(9) , the 
nonperturbative constant does not enter the calculation of LN model. [] But in non- Abelian case , 
especially including gluon self-interactions, as showed in Fig. 4 and Fig. 5, g connot be all absorbed by 
normalized condition Eq.(9), so g enters indeed into our the calculation. Unfortunately our knowledge 
for the non-perturbative coupling constant g is very poor now , in order to get a sensible answer for 
the gluon structure function at small x Q a n — j- of order one is needed, here we take the value 

Putting these parameters into EqjlT], the double diffractive production cross section of Eq.(l) for y/s 
from 20GeV to 2 x 10 4 GeV, where \t%\ <lGeV 2 , S — 0.1 are evaluated, as shown in Fig. 7. We see for 
yfs = 23.7GeV, 29Gev, 630GeV(SppS) and 1.8 TeV (Tevatron) energies, the double diffractive production 
cross sections are 1.6 x 10 2 n6, 2.5 x 10 2 nb, 2.8/i> and 4.6/A respectively. In the Joint CERN-IHEP experiment 
in 300GeV central ir~N collisions (\/s — 23.7GeV ),the production cross section in the range < xp < 0.3 
(the experiment setup lead to the acceptance is zero for xp < 0, ) is 0.2 ± 0.1/1& . Taking account of the 
addtivc quark rule, this corresponds to 0.4 — 1.4/i6 in P — P(P) central collisions at the same center-of-mass 
energy in the range —0.3 < xp < 0.3 . Since /(1500) is produced dominately at small \t\ Q][^5|, our results 
are reasonble . 

We have restricted ourseleves to calculate elastic diffractive production process. It is easily extend to 
diffractive dissociation processes by removing the form factor F\(ti)(i = 1,2) in Eq.(2) and thereafter , the 
resulting cross sections are also show in Fig. 7. We can see that as center of mass energies increased the ratio 
of the cross section of diffractive dissociation to that of elastic one is reduced from 6 at yfs = 23GeV to 3 at 

= 20Tev . 

Throughout the calculation we work in Feynman gauge. If one were to work in another gauge the function 
D(q 2 ) must be substantially modified to ensure that physical observables are gauge independent . 

In connection with the problem of experimentally detecting the process Eq.(l), we make following dis- 
cussions: 

1. Since glueballs in process Eq.(l) are produced through DPE in high energy p — p(p) collision, they 
should be laid predominately in the central region of rapidity distribution of final particles and with 
large symmetry rapidity gaps from each of the final proton ( antiproton ) direction. However several 
effects would be weaken the gap rate. First, althrogh the glueball and its decay product are colorless 
, as final state interacting there are get having survival soft color interactions between pre-glucball 
with outgoing pro-nucleons . The estimated gap survival rate to be 10~ 2 at Fermilab Tevatron p6| |. 
Secondly, the produced glueball is a light object (mo ~ 1.5GeV), they will be produced in a broad 
rapidity range, for example, at Tevatron those glueballs are in the rapidity region —5.9 < 77 < 5.9 . 

2. We take the D-L IP model and DPE process to discuss glueball production process at small \t\ , like 
all other high energy diffractive processed discussed by using the same model, the common remarkable 

3 Even in process which are absent for gluon selfinteractions, may meet the same troubles. For example, in discussiong the 
diffractive Higgs production process p + p-^p + p + H, only after supposing implicitely that coupling constant between top 
quark with IP equals to that of u, d quark with IP, then constant g can be all absorbed in Eq. (9). 
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character and parameter independent property are the energy dependence of total cross sections which 
are proportional to s 2 ( Q ^(°) _1 ) . So the cross section in this model always increases slowly with the 
increase of the center of mass energy, we could use this point to distinguish this production mechanism 
of glueball from others, especially from gluon - gluon fusion mechanism in parton model, since for later 
the production cross section would be decreased rapidly as energy increases, perhaps this is the most 
effective way to test our results if one could measure the total cross section of process Eq.(l) over a 
large energy range. 

In conclusion, using the field theroy model of pomeron exchange and the color singlet approximation of 
glueball, we obtain a parameter-free prediction of the cross section of ++ glueball in diffractive production , 
which combineing with the experiment measureing of a(hh — > hhG) ■ BR(G — ► hi + h2 + • • •) can output the 

important quantities BR(G — > hl + h2-\ ) , saying BR(G -> KK), or if we knew BR(G — > hi + h2 H ) 

then glueball production would be a test of the approach we used. 
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Figure Captions 



Fig.l. Sketch diagram for Eq.(l) via DPE. 

Fig. 2. G production in q — q(q) scattering via DPE, dotted lines are non-perturbative gluon of IP. 
Fig. 3. s— channel discontinuity of the first diagram in Fig. 2. q\, q2 lines arc limited on their mass shell. 
Fig. 4. Vertex of glueball and gluons of IP. Vertical dashed line limit cutting lines on mass shell. 
Fig. 5. The same as Fig. 4, but for four gluons vertex. 
Fig. 6. Glueball production in radiative decay of J/ip. 

Fig. 7. Production cross section of ++ glueball in double diffractive P — P(P) collision, a) Elastic case, b) 
Including the diffractive dissociation contribution 
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Fig.l Sketch diagram for Eq.(l) via DPE. 
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Fig. 2 G production in q 
perturbative gluon of IP. 



q(q) scattering via DPE, dotted lines are non- 
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Fig. 3 s— channel discontinuity of the hrst diagram in Fig. 2. qi } q 2 lines are 
limited on their mass shell. 



1 



h,fi,b 



a, e 




+ (a ^ (3; e ^ h) 



/ P,h . 
^ i 

.4 Vertex of glueball and gluons of IP. Vertical dashed line limit cutt 

lines on mass shell. 
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Fig. 5 The same as Fig. 4, but for four gluons vertex. 
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Fig. 6 Glueball production in radiative decay of J/ip. 
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